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1.  INTRODUCTION 


Airborne  particles  scatter  and  absorb  electromagnetic  (EM)  waves  and  affect  the  performance  of  target 
acquisition  systems  The  extinction  due  to  particles  depends  primarily  on  the  ratio  of  the  particle  size  and 
the  EM  wavelength.  Because  the  wavelengths  used  in  the  Army  target  acquisition  systems  and  other 
remote  sensors  range  from  visible  to  infrared  optical  wave  bands  and  millimeter  radar  waves  to  low- 
frequency  acoustic  waves,  the  particle  sizes  that  need  to  be  monitored  range  from  those  of  gas  and  aerosols 
to  coarse  sand.  This  research  project  seeks  to  develop  an  accurate  model  for  predicting  the  number 
concentration  and  size  distribution  of  airborne  sediment  particles  (including  sand  and  dust)  from  saltation 
and  suspension  caused  by  turbulent  wind. 

The  transport  problems  have  been  investigated  by  researchers  in  many  different  fields,  including 
atmosphenc  science,  geological  science,  agncultural  engineering,  chemical  engineering,  mechanical 
engineering,  ocean  engineering,  and  civil  engineering.  The  amount  of  research  results  and  literature 
references  is  extensive  and  diversified.  Much  of  the  work  has  been  reviewed  by  Bagnold  (1941),  Raudkivi 
(1976),  Yalin  (1977),  Clift  et  al  (1978),  Hwang  vl982),  SIcath  (1984),  Greeley  and  Ivcrscn  (1985), 
Anderson  (1986),  Werner  (1987,  1990),  and  Dyer  and  Soulsby  (1988),  among  others.  The  matenal 
presented  in  this  section  represents  merely  a  small  drop  in  a  very  big  bucket. 

In  general,  there  are  three  major  modes  of  sediment  transport:  creepuig,  saltation,  and  suspension.  For 
wind  (eolian)  transport  of  .sand  on  land,  saltation  is  of  special  unportance  because  it  dominates  a  rather 

thick  layer  (of  the  order  of  1  meter,  Bagnold,  1941)  of  the  lower  boundary  Suspension  is  important  m  the 

transport  of  fine  dust  and  aerosols,  which  can  be  airborne  for  a  long  period  of  time  (of  the  order  of  days  or 
weeks,  Toba,  1965)  Creeping  is  confined  in  the  immediate  vicinity  of  tJie  ground  and  contributes  to  a 
rather  small  fraction  of  the  total  transport  (Bagnold,  1941;  Nickling,  ,1978).  From  careful  observations, 
Bagnold  (p  10,  1941)  offered  the  followuig  description  of  salDtion  and  suspension  in  the  open  field: 

"A  popular  misconception  exists  regarding  sand  storms,  due  to  a  failure  to  distinguish  sand 
from  dust.  When,  in  any  arid  country,  after  a  spell  of  calm  weather,  a  strong  wind  begms  to 
blow  from  a  new  direction,  the  air  becomes  charged  with  a  mist  of  small  particles  WTere  the 
surface  is  alluvial,  with  little  or  no  sand  on  it,  such  as  in  Iraq  or  the  country  round  Khartoum, 

the  dust  rises  in  dense  clouds  to  a  height  of  several  thousand  feet  and  the  sun  is  obscured  for  a 

long  period.  This  is  obviously  a  dust  storm,  though  it  is  often  wrongly  designated  by  the 
possibly  more  thnlJing  and  cleanly  term  'sand  storm’.  Owing  to  their  feeble  terminal  velociiy 
of  fall  the  very  small  dust  particles  are  raised  and  kept  aloft  by  the  upward  currents  of  the 
wind's  internal  movements  ... 

On  the  other  hand,  in  an  erosion  desert,  the  only  free  dust  consists  of  those  fine  rock  particles 
which  have  been  loosened  by  weathering  since  the  last  wind  blew,  and  have  therefore  not  been 
earned  away.  In  such  country  the  wind  produces  for  the  first  hour  or  so  a  mist  consisting  of 
both  dust  and  sand.  Later,  although  the  wind  may  have  shown  no  signs  of  slackening,  the 
mist  disappears.  But  the  sand  still  continues  to  drive  across  country  as  a  thick,  low-flying 
cloud  with  a  clearly  marked  upper  surface.  The  air  above  the  sand  cloud  becomes  clear,  the 
sun  shines  again,  and  people's  heads  and  shoulders  can  often  be  seen  projecting  above  the 
cloud  as  from  the  water  of  a  swimimng-bath  Where  the  ground  is  composed  of  coarse  grains, 
pebbles,  or  large  stones,  the  top  of  the  cloud  may  be  2  metres  above  it,  but  is  usually  less. 

Where  the  surface  consists  of  fine  sand,  such  as  that  of  a  dune,  the  height  of  the  sand  cloud  is 
noticeably  lower" 


The  description  above  clearly  indicates  the  possible  senous  effects  of  eolian  processes  in  the  operation  of  a 
target  acquisition  system  and  other  remote  sensing  devices.  It  is  noted  that  the  physical  mechanisms 
governing  the  aerosol  transport  process  arc  very  complicated,  as  reflected  from  the  fact  that  the  sources  of 
aerosol  particles  can  be  from  many  thousand  miles  away,  there  is  no  reliable  model  to  desenbe  the  aerosol 
motion;  and  our  understanding  of  aerosol  distnbutions  relies  on  experimental  measurements.  Results  on 
the  size  distnbution  of  aerosols  have  been  generally  classified  mto  major  categories  such  as  Urban,  Rural, 
Remote  continental.  Desert  dust  storm.  Polar,  Stratosphere,  Sea  salt  and  Background  (e  g.,  Junge,  1972, 
Jaenicke,  1988).  For  larger  particles,  including  sand  grains  and  dust  particles,  the  distribution  is  more 
localized  and  amenable  (in  a  relative  sense)  to  mathematical  treatment.  The  capability  to  establish  the 
transport  rate  or  concentration  profiles  as  functions  of  wind  and  particle  parameters  underscores  the 
significant  advance  of  our  understanding  of  the  large  particle  dvriamics  (although  admittedly  still  far  away 
from  perfect  or  even  satisfactory). 

W’ltli  the  advance  of  digital  computers,  many  numerical  models  have  been  developed.  Mechanisms  such  as 
sediment-wmd  mtcraction  and  the  random  nature  of  parameters,  including  particle  inception  by  wind  force, 
variabiliiy  of  ejection  angles,  impact  angles,  number  of  particles  mobilized  by  unpact,  and  panicle 
trajectories  subjected  to  drag,  lift,  and  spinning  effects,  can  be  incorporated  (c  g  .  Owen,  1964,  White  and 
Schulz,  1977,  .Anderson,  1986,  .Anderson  and  Hallet,  1986,  Ungar  and  Haff,  1987,  .Anderson  and  Haff, 
1988,  1990,  Haff  and  Anderson,  1989,  Werner,  1990).  These  models  enhance  our  understanding  of  the 
physics  of  sediment  transport  considerably.  It  is  noted  that  such  numencal  models  are  computationally 
intensive  and  may  take  hours  or  even  days  to  complete  a  "run"  of  relatively  sunple  scenario  (Haff,  personal 
communication,  1991);  thus,  the  applicability'  to  field  use  remams  luruted. 

Although  practical  application  using  results  from  modem  numerical  sunulations  on  the  physics  of  sedunent 
transport  appears  to  be  remote,  there  are  some  robust,  time-tested  relationships  denved  from  sediment 
research  that  deserve  a  closer  examination.  One  of  the  most  significant  developments  is  the  establishment 
of  sediment  transpon  rate  equations,  formulated  to  predict  the  flux  of  sediment  as  a  function  of  wind  and 
particle  properties.  The  development  of  these  formulas  was  based  on  the  consideration  of  momentum 
balance  between  airborne  sedunent  and  wind  motion,  as  vsell  as  data  fitting  to  account  for  physical 
mechanisms  that  were  not  easily  formulated  mathematically  (Bagnold,  1941).  Since  Bagnold’s  rmlestone 
work,  numerous  new  equations  have  been  mtroduced,  but  tune  and  again  the  fundamental  physical 
mechanism  described  by  Bagnold  was  found  to  be  rather  accurate. 

The  transport  rate  equation,  however,  provndes  the  total  sediment  flux  only  and  offers  very  poor  size 
resolution.  As  mentioned  earlier,  for  remote  sensing  applications,  a  good  understanding  of  the  size 
distnbution  of  airborne  particles  is  cniical  because  the  properties  of  absorption,  scattenng,  and  ex-tmction 
are  wavelength-  and  particle-size-dependent  (e  g.,  Mie,  1908,  Rayleigh,  1945,  Sonunerfcid,  1959;  van  de 
Hulst,  1981).  In  order  to  apply  the  results  of  sedunent  transport  research  to  target  acquisition  applications, 
we  need  to  develop  a  sequence  of  mathematical  manipulations  to  bnng  out  the  mherent  information  of  the 
size  distribution  function  of  airborne  particles.  This  is  done  by  introduemg  an  equivalent  concentration  to 
represent  the  saltation  particles.  By  equating  the  transport  rate  with  the  rate  of  sediment  flux  (the  product 
of  concentration  and  convective  velocity)  over  the  saltation  layer,  it  is  possible  to  recover  the  mformation  of 
the  size  distribution.  W'ith  the  size  distribution  known,  the  number  of  airborne  particles  within  the  path  of  a 
target  acquisition  system  can  be  calculated. 

For  smaller  particles  that  may  suspend  in  the  air  once  lifted  from  the  ground,  the  dominant  process  is 
turbulent  diffusion.  With  proper  assumptions  of  eddy  diffusivity,  the  resulting  concentration  profiles  decay 
with  distance  from  the  bottom  bound  iry  foliov^^^g  either  a  power  law,  an  exponential  distribution,  or  some 
combined  form  of  both  (eg.,  Hwang,  1982).  The  suspension  and  saltation  processes  arc  closely  related  In 
particular,  the  impact  of  saltation  particles  causes  the  dislodging  of  small  dust  particles  to  start  the 
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suspension  process,  and  the  saltation  layer  controls  the  boundary  condition  of  the  suspension  layer.  Tlie 
number  densiy  and  size  distribution  of  the  saltation  layer  determine  the  absolute  concentration  of  the  much 
thicker  suspension  layer.  Coupling  of  these  tv,'o  layers  is  important  in  denying  the  total  density  of  particles 
under  any  given  configuration  of  a  target  acquisition  system. 

In  order  to  interface  with  the  operations  of  target  acquisition  systems,  a  procedure  for  coordirute 
conversion  from  Cartesian  coordinates  to  spherical  coordinates  is  developed.  Computational  examples 
show  that  the  integrated  number  of  particles  within  a  pencil  beam  of  a  remote  sensor  increases  significantly 
as  the  incident  angle  increases;  a  few  orders  of  magnitude  increase  from  10®  to  80®  is  common  for  a  typical 
vertical  stratification  of  airborne  particle  distribution. 

This  report  documents  a  Phase  I  study  to  evaluate  and  develop  models  of  the  saltation  and  suspension  of 
sediment  by  turbulent  transport.  In  Section  2,  a  brief  discussion  is  presented  of  the  theoretical  background 
of  saltation,  suspension,  the  flux  formulation  to  retrieve  the  size  distribution  information,  and  the 
mathematical  aspect  of  coordinate  transformation  from  Cartesian  to  sphencal  coordinates  In  Section  3,  a 
comparison  and  evaluation  is  given  of  a  small  group  of  sediment  transport  rate  equations  The 
parameterization  scheme  to  facilitate  numencal  evaluation  of  number  density  from  the  flu.x  formulation  is 
given  in  Section  4.  Results  and  discussions  of  sensitivity  tests  performed  to  investigate  and  rank  the 
importance  of  numerous  parameters  are  given  m  Sections  5  and  6,  respectively,  and  recommendations  for 
further  research  are  presented  in  Section  7. 
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2.  THEORETICAL  BACKGROUND 

Two  modes  of  sediment  transport  are  of  special  interest  in  the  study  of  airborne  particle  distribution, 
suspension  and  saltation  The  former  governs  the  dispersion  of  smaller  dust  particles  into  the  air,  and  the 
latter  governs  the  ejection  of  larger  sand  grains  from  the  surface.  The  impact  of  sand  grains  is  an 
important  mechanism  for  dislodging  the  smaller  dust  particles  that  are  either  shadowed  by  larger  grains  or 
present  in  the  source  material  as  aggregates. 


2.1  Saltation 

One  of  the  most  significant  relationships  derived  from  sediment  transport  research  is  the  establishment  of 
the  transport  rate  as  a  function  of  the  wind  friction  velocity  or  wind  speed.  Because  the  contribution  of 
different  size  classes  to  the  total  sediment  transport  is  proportional  to  the  volumes  of  Jie  size  classes,  the 
transport  rate  is  dominated  by  the  larger  particles,  which  move  in  the  air  in  saltation.  For  e.xample, 
Bagnold  (194 1)  reported  that,  for  regular  desert  sand,  the  saltation  load  is  80%  of  the  toDl  sediment  flu.\. 
To  establish  a  relationship  be^veen  the  transport  rate  and  the  turbulent  vsind  condition,  Bagnold  (1941) 
considered  the  momentum  balance  of  airborne  particles.  Startmg  with  a  saltating  particle  of  mass  m  that 
nses  from  the  surface  wnth  a  honzonDl  velocity  u,  and  strikes  on  the  surface  agam  after  traveling  a 
distance  /  with  a  horizontal  velocity  and  assuming  that  all  of  the  honzontal  momentum  was  lost  after 
impact,  the  averaged  momentum  extracted  from  the  air  per  unit  length  is  therefore  Extendmg 

the  concept  of  one  grain  to  the  mass  flux  of  sand  q,  the  rate  of  momentum  loss  of  air  becomes  q{u2~u^),'l. 
This  momentum  loss  results  in  a  resisting  force  r  on  the  surface,  if  u,«ii2,  and  the  following  expression  is 
established; 


q^=T  =  puJ^  (1) 

where  p  is  the  density  of  air  and  u,  is  the  wind  friction  velocity.  From  grain  path  analysis,  Bagnold  further 
found  that  over  a  wade  range  of  values,  where  w,  is  the  initial  vertical  velocity  of  rise  of  the 

gram  at  the  beginning  of  its  path  and  g  is  the  gravitational  acceleration.  The  sediment  flux  becomes 

q  =  (2) 

g 


Assuming  that  w,  is  proportional  to  u,,  based  on  the  argument  that  the  grains  are  exposed  to  the  air  stream 
at  all  heights  and  that  their  mean  velocity  of  unpact  is  likely  to  be  controlled  by  u,  which  defined  the  wind 
speed  distribution,  Bagnold  reached  the  conclusion  that  the  rate  of  sediment  flux  q  is  proportional  to  the 
cubic  power  of  the  wind  friction  velocity; 


g 


(3) 


where  a  is  the  representative  grain  diameter  of  the  source  sediment,  =  0.025  cm  is  a  reference  sand  size 
and  fl  is  an  empirical  coefficient  determined  to  be  1.5  for  a  nearly  uniform  sand,  1.8  for  naturally  graded 
sands,  such  as  those  found  on  sand  dunes,  and  2.8  for  a  sand  with  a  very  wide  range  of  grain  sizes.  A  great 
number  of  newer  equations  have  been  suggested  over  the  past  few  decades  (e  g.,  see  Greeley  and  Iveisen. 
1985,  Sarre,  1987),  but  the  fundamental  physics  of  momentum  balance  between  sand  flow  and  wind 
motion  has  not  changed,  and  nearly  all  of  the  formulas  maintain  the  cubic  power  law  relationship  (either 
directly  with  wind  friction  velocity  or  the  velocity  scale  representmg  the  excess  wind  stress,  which  is  the 
difference  of  wind  stress  and  particle  threshold  stress). 
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2.2  Suspension 

Particle  suspension  is  governed  by  turbulent  diffusion  in  the  fluid.  This  is  generally  formulated  as  the 
balance  of  turbulence  agitation  with  flow  convection  and  gravitational  settling.  Efforts  to  study  the 
suspension  problem  have  been  performed  by  many  scientists  in  various  fields,  such  as  atmospheric  science, 
hydraulics,  oceanography,  coastal  engineering,  chemical  engineering,  agriculture,  and  environmental 
science.  The  process  can  be  expressed  as  (e  g.,  Yalin,  1977;  Hwang,  1982;  Anderson,  1986) 

“  +  V«(u.c)  =  V«(vVc)  (4) 

where  c  is  the  (volumetric)  concentration  of  the  suspended  particles,  Up  the  convective  velocity  vector  of  the 
panicles,  vthe  fluid  (air)  viscosity  and  r  the  time.  It  is  rximmon  to  express  the  variables  c  and  u  in  two 
components:  the  mean  and  the  fluctuation  (due  to  the  turoulent  nature  of  the  process),  that  is,  c=C'+c'  and 
Up  =  Up  +Up  .  The  equation  governing  the  mean  concentration  becomes 

^  +  V.((CUp)  +  (c'up'))  =  V.(vVC)  (5) 

where  the  angle  brackets  represent  the  mean  value  (ensemble  or  temporal  average)  of  the  bracketed  prop- 
eny.  A  turbulent  diffusivity  coefficient  is  usually  introduced  to  represent  the  diffusion  term  of  the 
turbulent  "flux": 

„  /  A  X  X 

(6) 

where  and  are  the  turbulent  diffusivitN'  in  the  x*.  y-,' and  z-directions.  The  solution  of  Eq  (5) 

depends  ver>'  much  on  the  turbulence  properties  of  the  flow  field.  Due  to  our  limited  understanding  of  flow 
turbulence,  closed-form  solutions  are  available  only  for  a  few  simplified  conditions.  For  example,  under 
the  assumptions  of  temporal  steadiness  and  honzontal  homogencitx  ,  also  neglecting  the  term  of  molecular 
diffusion  ui  companson  with  the  turbulent  diffusion  terms,  the  resulting  equation  becomes 

where  wyis  the  fall  velocity  of  the  particle.  The  solution  of  Eq  (7)  is 


C(r)  =  Qexp  -[  -^dz 


in  which  is  the  reference  concentration  at  elevation  Zg.  Depending  on  the  definition  of  the  turbulent 
diffusivity',  the  distribution  follows  either  an  exponential  (for  consiam)  or  power-law  (for  proportional 
to  z)  decay  with  distance  froni  the  boundary.  Other  more  complicated  definitions  of  would  lead  to  more 
complicated  vertical  profiles  of  suspension.  The  rapid  decay  from  the  boundary  remains  the  common 
feature  and  stresses  the  importance  of  including  the  vertical  variation  of  dust  and  aerosol  distribution  when 
the  results  are  to  be  applied  close  to  the  boundary. 

When  the  bottom  material  is  not  composed  of  uniform-size  particles,  the  size  distnbution  of  the  source 
matenal  needs  to  be  taken  into  consideration.  If  the  interaction  bertveen  different  size  classes  is  negligible, 
linear  superposition  can  be  applied  and  the  resulting  total  concentration  is 


C'(z)  =  j  p{a)C(a,z)da 


5 


where  p(a)  is  the  size  distribution  of  bottom  material  and  C(a,  z)  is  the  suspension  of  an  individual  size 
class,  which  takes  the  fonn  of  Eq.  (8)  with  the  fall  velocity  and  turbulent  diffusivity  properly  determined 
for  each  size  class. 


2.3  Sediment  Flux  and  Size  Distribution 


From  empirical  obser/ations,  it  is  generally  found  that  the  ratio  of  saltation  and  suspension  transport  is 
approximately  constant.  Bagnold  (1941)  found  that  the  ratio  is  4:1  based  on  laboratory  measurements 
using  sand  ranging  in  size  from  0. 1  to  1 .0  mm.  Nickling  (1983)  conducted  field  measurements  on  the  Slims 
River  Valley,  Yukon  Territory,  where  the  mean  sediment  size  is  considerably  finer  (0.06  to  0.12  mm)  and 
the  fraction  of  suspension  increases  to  as  much  as  69%  of  the  total  transport.  Regardless  of  the  exact 
proportionality  of  saltation  versus  suspension,  the  sediment  transport  rate  can  be  expressed  in  terms  of  the 
airborne  sediment  concentration,  C(a,  z),  here  representing  the  sum  of  suspension  and  saltation.  The 
volume  concentration,  C(a.  zj.  is  further  related  to  the  size  distribution  of  the  suspension,  n{a.  z),  by 


C(a,z)  =  nia,z) 
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(10) 


By  definition,  the  transport  rate  can  be  wrinen  as  the  vertical  integration  of  the  product  of  sediment 
concentration  and  convective  velocity,  i.e.. 


<7  = 


(a,z)  da  dz 


(11) 


Combinii.g  Eqs.  (3)  and  (11),  it  is  immediately  obvious  that  the  size  distribution  as  a  function  of  wind 
condition  and  iinface  material  property  can  be  obtained  from  the  sediment  transport  formula: 


4/ta^  ( a 

-^Up{a,2)dadz=  b\ 


(12) 


when  the  integration  limits  are  explicitly  expressed  as  (o,,  Oj)  and  Zj),  respectively.  Equation  (12) 
contains  many  unknowns  including  the  size  distribution  function  n{a,  z),  the  convective  velocity  i/^,  and  the 
four  vanables  defining  the  uitcgration  limits.  Our  pnmary  tasks  are  to  provide  additional  relationships 
based  on  theoretical,  numerical,  and  expeiomental  results  so  that  enough  equations  can  be  formulated  to 
solve  for  the  size  distnbution  function.  For  example,  the  integration  limits  (o,,  Oj]  represent  the  range  of 
panicle  sizes  that  can  be  mobilized  under  a  given  wuid  condition,  which  can  be  denved  from  our  knowledge 
of  the  threshold  criterion  of  particle  motion,  and  the  factor  Zj  defines  the  upper  limit  of  uitegration,  which  is 
properly  interpreted  as  the  saltation  height.  More  discussion  of  the  parameterization  is  given  in  Section  4, 


2.4  Spherical  Coordinates  Transformation 

The  sediment  concentration  or  mass  flux  formulas  are  generally  developed  in  tf<e  rectangular  coordinate 
system.  To  predict  the  concentration  of  airborne  sand  and  dust  along  the  path  between  observer  and  target, 
it  is  more  useful  to  express  these  terms  in  spherical  coordinates.  Fig,  1  shows  the  two  coordinate  systems 
from  which  the  following  relationships  can  be  established: 


r  =  sin  cos  0 

(13) 

>'  =  /?  sin  sin  t? 

(14) 

z  =  p  cos^  +  D 

(15) 
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Figure  1.  Sketch  depicting  the  wave  path  of  a  target  acquisition  system  in  rectan¬ 
gular  Cartesian  coordinate  system  and  spherical  coordinate  system. 

(Note  tlie  svmbols  used  here  follow  the  nruthematical  convention  with  ^  as  the  zenith  angle  and  ^  as  the 
azimuthal  angle.  In  radar  applications,  6  and  ^  are  usually  reversed,  that  is,  9  is  the  incident  or  zenith 
angle  and  ^  die  azimi.thal  angle).  The  vertical  of.set  D  is  the  elevation  of  the  target  acquisition  system 
The  integration  region  of  interest  is  the  path  cone  s.aown  in  the  figure.  Because  the  extinction  parameter  of 
the  urget  acquisition  system  is  sensitive  to  the  instrument  wavelength  and  particle  size,  the  desired  quantity 
from  the  solution  algorithm  is  the  concentration  of  particles  of  a  given  size,  that  is,  the  size  distnbution 
n{a;  p.  6.  4>)  at  different  ranges  and  azimuth  and  zenith  angles.  With  the  size  distribution  given  in  the 
rectangular  coordinates,  the  transformation  to  spherical  coordinates,  and  the  integration  to  obtain  the  inte¬ 
grated  size  distribution  within  the  path  cone  can  be  obtained  using  Eqs.  (13)  through  (15).  As  discussed 
earlier,  we  have  assumed  horizontal  homogeneity  in  the  airborne  particle  distribution;  thus,  the  only 
fiinction  involved  in  the  coordirute  transformation  is  the  vertical  distribution  function  Following  the  nou- 
tion  to  be  used  in  Section  4,  we  can  wnte  n{a;  p,  0,  ^  =  n(a,  z)  =  Kn{a)J[a,  z),  where  y(fl;  z)  =j{a,  pzoitp) 
is  the  vertical  distribution  fiinctic The  number  of  particles  integ  tted  over  the  path  cone  is  therefore 
proportional  to 

F(a)  =  JJJ/(o;z)  dxdydz  =  jjj/ia  ,pcosf>  s\n^dpd(j)dO 

MOI  HD/ 


where  p^sin  ^  is  the  Jacobian  of  the  coordinate  transformation,  F[a)  is  the  integrated  function,  and  vol  de¬ 
notes  the  volume  of  the  air  space  within  the  wave  path,  as  defined  in  Fig.  1: 


SO  SO  d)n  ~  S9  in 
vol  =\Q< p<R,  — 2 — 


(17) 


Using  the  solid  angle  expression  as  in  remote  sensing  applications,  Eq  (16)  can  be  written  as 

F(a)  =  jjj/ia, pcosS)/^  dCl  (jg) 

¥CJ 

where  dfi=  sin  ^  c/d  is  element  of  the  solid  angle.  The  functional  dependence  of/(d)  suggests  that 
the  inclusion  of  airborne  particles  in  the  remote  sensing  computation  is  equivalent  to  introducing  a  modifi¬ 
cation  of  the  antenna  pattern  (e.g.,  Kerr,  1951;  Stutzman  and  Thiele,  1981).  The  transformation  of  the 
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cross-sectional  area  and  volume  concentration  of  particles  in  the  path  cone  can  be  derived  in  a  similar 
fashion  as  that  of  the  number  density,  Eq.  (16). 

For  a  smooth  function /such  as  generally  encountered  in  the  vertical  distribution  of  airborne  particles,  the 
integral  can  be  carried  out  and  a  closed-form  solution  is  available.  For  example,  if  the  vertical  distributic  i 
of  airborne  particles  follows  a  power-law  decay  form,y(/)  =  r”*,  the  integrated  number  of  particles  within  a 
pencil  beam  of  a  target  acquisition  system  is  shown  in  Fig.  2  as  a  function  of  the  incident  angle.  In  this 
figure,  the  computed  number  is  normalized  such  that  when  the  particle  distribution  is  uniform  (m  -  0),  the 
integrated  number  is  I .  As  clearly  illustrated  in  the  Egure,  at  a  low  grazing  angle  (high  incident  angle),  the 
encountered  number  of  airborne  particles  can  increase  drastically,  more  than  lOOO-fold  between  10°  and 
85°  for  m  ^  2.9.  This  large  increase  of  particle  concentration  at  the  grazing  angle  can  significantly  affect 
the  performance  of  target  acquisition  systems  in  near-ground-level  applications. 

2.5  Research  Tasks 

Four  research  tasks  were  identified  for  the  Phase  I  effort: 

•  Task  1.  Evaluation  and  selection  of  sediment  transport  formulas. 

•  Task  2.  Derivation  of  size  distribution  from  sediment  flux  formulation 

•  Task  3  Establishment  of  parameterization  scheme  to  facilitate  numerical  computation 

•  Task  4  Performance  of  sensitivity  tests  and  evaluation  of  the  feasibility  of  the  proposed  concept  to 

derive  the  size  distribution  function  from  sediment  transpon  rate  equations. 


Incident  Angle 

Figure  2.  The  effect  vertical  distribution  of  airborne  particles  on  variation  of  the 
Integrated  particle  concentration  within  the  pencil  beam  of  a  target  acqui¬ 
sition  systf  m  plotted  as  a  function  of  Incident  angle.  The  vertical  decay  func¬ 
tion/?)  is  assumed  a  power-law  of  the  form  r".  The  number  density  for  a  unifonn 
particle  distribution  (m=0)  is  used  as  the  normalization  factor.  Tlie  normalized  curves 
for  m=0,  0  5,  I.l,  2,  and  2.9  arc  shown. 
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3.  COMPARISON  OF  SEDIMENT  TRANSPORT  RATE  EQUATIONS 

Seven  sediment  transport  rate  formulas  (Eqs  (19)  to  (25)  given  in  Table  1]  are  compared  to  evaluate  the 
statistical  difference  among  these  equations.  The  common  feature  of  these  equations  is  the  cubic-power 
wind  dependence,  either  in  terms  of  the  wind  friction  velocity,  u.,  or  the  wind  speed  at  a  given  reference 
level,  f/,00  (taken  as  1  meter  above  the  ground  in  this  report).  The  more  recent  equations  usually 
incorporate  the  threshold  friction  velocity  in  the  formulation  to  emphasize  that  transport  occurs  only  when 
the  surface  stress  exceeds  the  resistance  of  grains  to  wine  forcing.  Formulas  of  Bagnold  (1941)  and  Zingg 
(1952)  do  not  contain  although  the  corresponding  expressions  using  wind  speed  U  have  a  threshold 
wind  speed  factored  in  Bagnold's  case  (p,  69,  Bagnold,  1941).  These  equations  are  understood  to  be 
applied  for  u,>u^.  In  this  report,  numerical  computation  based  on  the  seven  transport  rate  equations  is 
performed  for  wind  speeds  ranging  from  2  to  40  m/s  and  characteristic  sediment  diameters  from  0.005  to 
0.05  cm  (50  to  500  /m\).  The  results  are  presented  in  this  section. 

Table  1.  Selected  sediment  transport  rate  equations  evaluated  in  this  study. 


B.Al 

Bagnold  (1941) 

/  A®* 

^  =  18  —1 

(19) 

KA 

Kawamura  (1951) 

q  =  i.n 

/ 

1  + 

\ 

m 

(20) 

l&L 

Lettau  and  Lettau  (197$) 

q  =  2.85 

1  j 

P,  3 

—  U. 

g 

(21) 

IS 

iversen  et  al.  (1976), 

Schmidt  (1982) 

<y.0.4^fl-^| 

u.,  1,  u.  ) 

J 

~u, 

8 

(22) 

OW 

Owen  (1964) 

mBM 

(23) 

21 

Zingg  (1952) 

^  =  0.83 

^^0  J 

0.75 

g 

(24) 

WH 

White  (1979) 

m 

1-^ 

{  «•  )g 

(25) 

3.1  Fundamental  Parameters 

Most  of  the  transport  rate  equations  are  expressed  in  terms  of  the  wind  friction  velocity,  u,,  and  the 
sediment  diameter,  a.  The  latter  is  one  of  the  most  important  parameters  of  the  bottom  material.  If  the 
material  is  uniform  (monodisperse),  the  dcfmition  of  sediment  size  is  straightforward.  For  mixed  materials 
(polydispersc),  the  median  size  is  usually  used  Two  other  parameters  that  are  also  included  in  some  of  the 
transport  rale  equations  arc  the  threshold  velocity,  iv.,,  and  the  particle  settling  velocity,  The  functional 
expressions  of  u.,  u.,,  and  nyarc  discussed  in  the  following. 
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Wind  Friction  Velocity 

Tlie  vsind  fnction  velocir>’  is  a  measure  of  ihe  wind  stress  applied  on  the  sediment.  If  the  bottom  material  is 
not  moving,  u,  and  U(z)  are  connected  by  the  well-established  logarithmic  distribution  (Bagnold,  1941; 
Schlichting,  1968) 

U(z)  =  5.15iu\ogj  (26) 

where  it  is  a  length  scale  characterizing  the  surface  roughness.  For  a  flat  surface  covered  with  uniform 
sand,  k  =  a/30.  The  wind  friction  velocity  can  be  calculated  from  the  reference  wind  velocity  at  1  meter, 

^100* 


751ogiHKlM  ”  19.99 -5.75 log  a 


The  vertical  distribution  of  wind  is  modified  by  sand  motion,  as  investigated  in  detail  by  Bagnold  (1941) 
The  resulting  velocity  profile  is  given  as 

U{:)  =  5.15u»log—*U,  (28) 

k‘ 

where  the  modified  roughness,  k,  and  threshold  velocity,  are  functions  of  particle  size  The  following 
formulas  (in  cgs  units)  are  adopted  following  the  laboratory  results  of  Zingg  (1952;  see  also  Honkawa, 
1978;  Horikawa  et  al.,  1986): 

lOrr  (29) 

U,  =  8940a  (30) 

The  wand  friction  velocity  on  a  movable  boundary  becomes 

_  ^100 -89400  _  t/,oo-8940o 

"  5  751orM  "5.75(1 -logo) 

^  lOo 

Again,  the  wind  velocity  at  1  meter  is  used  as  a  reference.  The  wind  friction  velocity  with  a  moving  sedi¬ 
ment  bottom  is  considerably  larger  than  that  on  a  fixed  rough  boundary  Fig  3  displays  the  computed 

results  ofEqs.  (27)  and  (31)  for  sediment  size  o=0.005,  0.01 .  0  05  cm  The  fnction  velocity  on  moving 

sand  (dotted  curves)  is  nearly  twice  the  magnitude  of  the  fixed  boundary  results  (solid  curves)  at  higher 
wind  velocities  (f/ioo^iO  m/s). 

Bagnold  also  conducted  field  measurements  on  natural  sand  dunes  and  found  that  for  standard  sand  dunes, 
^'*=  1  cm  and  (7,  =  4  m/s.  With  these  values,  the  wind  fnction  velocity  would  be 

u.  =  0.087(7, 00 -34.7  (32) 

The  calculated  results  using  Eq.  (32)  are  shown  as  x's  in  Fig.  3  also.  It  appears  that  Eq  (31)  derived  from 
laboratory  results  may  underestimate  the  friction  velocity;  however,  lacking  a  more  comprehensive  data  set 
from  field  measurements,  it  is  used  in  the  .ollowing  evaluation  of  various  sediment  transport  equations. 


Threshold  Friction  Velocity 

As  wind  blows  on  a  sandy  surface,  it  exerts  drag  on  sand  grains.  Due  to  gravitational  stabilization,  the 
sand  gram  resists  the  drag  until  wind  exceeds  some  cntical  velocity  (the  thre.shold  velocity)  to  set  the  sand 
panicle  in  motion.  The  threshold  velocity  can  be  calculated  by  considering  tiic  balance  of  moments  due  to 


(cm/s) 


^100 

Figure  3.  Wind  friction  velocity  as  a  function  of  wind  velocity  at  1  •meter  elevation. 

Solid  curves:  rough  turbulent  boundary  layer  v.ith  fixed  roughness  elements  (uniform 
sand  of  diameter  from  0.005,  0.01, 0.05  cm),  u, increases  vxith  increasingo  for  the 
same  17, qq.  Dotted  curves:  movable  boundary  layer  with  uruform  sand  of  the  same 
diameter  range  as  previous  case,  x:  results  corresponding  to  the  field  observations  on 
standard  sand  dunes  by  Bagnold  (1941). 

the  two  competing  forces,  the  wind  drag,  Fp,  and  the  gravitational  force,  Fq  (Bagnold,  1941).  With 
reference  to  Fig  4a,  the  rw  o  forces  can  be  expressed  as 

Fo=^ifi-p,)g  (33) 

6 

(34) 

2  4 

where  A'  is  z  factor  representing  the  fraction  of  the  particle  cross>sectional  area  exposed  to  the  flow,  and 
Cp  is  the  drag  coefficient.  Alternatively,  u,‘^-CplF  can  be  substituted  into  Eq.  (34)  and  the  drag  force 
expressed  in  terms  of  the  friction  velocity.  The  corresponding  moments  about  the  resting  point  of  the  grain 
from  the  two  force  components  are 

A/c  =—(/?- sin  o  (35) 

6  Z 

Mn  =  —ai,^-^^—A'*—cosa  (36) 

2  4  2 


II 


Figure  4.  (a)  Sketch  showing  the  inception  of  particle  motion  subject  to  fluid  dis¬ 
turbance;  (b)  experimental  curves  of  the  threshold  friction  velocity  (from 
Bagnold,  1941,  Fig.  28);  (c)  sketch  showing  the  force  balance  of  a  particle 
falling  in  still  air;  and  (d)  settling  velocity  curve  (in  dimensionless  form, 
from  Yalin,  1977,  Fig.  3.2). 
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where  a  is  approximately  the  angle  of  repose.  At  the  critical  condition  Mq=Mq,  by  equating  (35)  with 
(36),  the  llireshold  fnction  velocity  is  found  to  be 

18  tana  Ip,  - p  .  \p,~ p 

The  experimental  value  for  A  in  the  case  of  sand  transport  in  air  is  found  to  be  approximately  0.1  for 
particles  larger  than  0.1  mm  in  diameter  (Bagnold,  1941;  Chepil,  1945)  For  smaller  particles,  the  grains 
are  not  exposed  to  the  wind  flow  directly  but  shielded  under  the  viscous  layer,  and  most  of  the  wind  drag  is 
absorbed  by  the  surface  rather  than  the  particles.  Subsequently,  the  threshold  friction  velocity  increases. 
Experimental  evidence  indicates  that  the  critical  diameter  of  such  transition  is  at  0.08  mm  (Bagnold,  1941; 
Chepil,  1945).  In  the  following  computation,  the  threshold  velocity  as  a  function  of  sediment  diameter  is 
digitized  from  Fig.  4b  (reproduced  from  Bagnold,  1941,  Fig.  28)  and  fitted  with  a  poKnomial  function  for 
digital  computation. 


Particle  Fall  Velocity 


Tlie  terminal  velocity  of  a  particle  falling  in  still  air  can  be  calculated  by  considering  the  balance  of  the 
grain  weight  and  the  drag  force  on  the  particle  (Fig  4c).  The  expressions  for  the  two  force  components  are 
the  same  as  those  given  m  Eqs.  (33)  and  (34)  except  that  ui  this  case  the  flow  velocity  is  the  fall  velocity, 
Wy;  and  the  panicle  is  fully  exposed  to  air  flow,  thus  the  factor  A'  is  unity  The  resulting  expression  for  the 
fall  velocity  is  _ 


-  M 

■V3C,,  p 


ga 


(38) 


The  drag  coefficient  is  a  function  of  the  Reynolds  number,  Re=Wja/v,  where  vis  the  fluid  viscosity.  In 
the  laminar  flow  regime  (/?e<l),  Cp  is  approximately  lAIRe,  and  wyis  proportional  to  oV  In  the  turbulent 
flow  regime  (/?e>1000),  the  drag  coefficient  is  approximately  constant,  and  the  fall  velocity  is  proportional 
to  the  square  root  of  the  gram  diameter  The  most  frequently  encountered  conditions,  howev  er,  fall  between 
the  two  regimes,  and  the  expression  of  Cq  is  rather  complicated.  Compiled  from  many  data  sets,  Yalin 
(1977)  constructed  a  dimensionless  plot  relating  the  fall  velocity,  grain  diameter,  and  fluid  viscosity  (Fig. 
4d).  The  curve  is  digitized  and  fitted  with  a  polynomial  function  for  numencal  computation. 


3.2  Results 

Seven  transport  rate  equations  [Table  1,  Eqs.  (19)  to  (25)J  were  compared.  The  computed  range  of  wind 
speed  at  1  meter  high  (C/jo«)  was  from  2  to  40  m/s,  the  wind  friction  velocity  (u.)  was  up  to  220  cm's 
(roughness  dependent),  and  the  range  of  sediment  diameters  (a)  was  from  0.005  to  0.05  cm  (50  to  500  /rm) 
The  complete  set  of  computed  results  is  presented  in  Figs.  5a  to  5j  and  6a  to  6g  The  following  conclusions 
were  reached  from  the  computed  results: 

•  Bagnold's  (1941)  and  Zingg's  (1952)  equations  do  not  contain  the  threshold  parameter.  Tliese  two 
equations  are  understood  to  be  applicable  for  winds  exceeding  the  threshold  condition  and  may  be  less 
accurate  at  lower  wind  velocities. 

•  In  general,  the  results  from  Bagnold's  formula  are  near  the  average  of  all  predictions  The  fonnulas  of 
Kawamura  (1951),  Lettau  and  Lettau  (1978),  and  White  (1979)  predicted  higher  transport  rates  than 
Bagnold's,  while  Zingg  (1952),  Owen  (1964),  Iversen  et  al.  (1976),  and  Schmidt  (1982)  predicted 
lower  rates. 
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(a)  10» 


(b)  lO-i 


(c)  10-1 


Figure  5.  Transport  rate  as  a  function  of  wind  friction  velocity  (left)  or  wind  velo¬ 
city  (right)  at  1 -meter  height.  Equations  (Tabic  1)  suggested  by  Bagnold  (solid 
curve);  Lettau  and  Lettau  (dotted  curve),  Kawamura  (dashed  curve),  Iversen  et  al.  and 
Schmidt  (circles);  and  Owen  (crosses)  arc  shown  for  uruform  sand  of  diameter  (in  cm) 
(a)  0.005,  (b)  0.01,  (c)  0.015,  (d)  0.02,  (e)  0  025,  (0  0.03.  (g)  0.035,  (h)  0.04, 

(i)  0.045,  and  (i)  0.05. 


Tiirrrrr 


•  The  computed  results  agree  with  each  other  within  a  factor  of  2  in  most  of  the  cases,  except  for  the 
formula  of  Owen  (1964),  which  may  be  off  by  a  factor  between  3  to  4  from  the  average  at  higher  wind 
velocities.  This  degree  of  agreement  is  considered  very  good  given  that  most  of  the  sediment  transport 
rate  data  show  order-of-magnitude  scatter  within  the  individual  data  set  and  even  more  when  different 
sets  are  combined  The  large  scatter  occurs  in  both  laboratory'  and  field  conditiof^s. 

•  The  expressions  of  Kawamura  (1951)  and  Lettau  and  Lettau  (1978)  produce  nearly  identical  results 
over  the  whole  range  of  conditions  tested.  The  former  is  somewhat  more  complicated. 

•  The  degree  of  size  dependence  varies  among  the  compared  equations  and  can  be  roughly  divided  into 
two  groups.  Four  of  the  seven  equations  show  distinct  size  dependence  (Fig.  6);  these  correspond  to 
the  equations  that  specified  the  particle  diameter  explicitly,  either  as  a  itself  in  the  cases  of  Bagnold 
(1941)  and  Zingg  (1952)  or  as  the  fall  velocity  iit  the  cases  of  Owen  (1964)  and  Iversen  et  al. 
(1976). 

•  In  the  other  three  formulas  (Kawamura,  1951,  Lettau  and  Lettau,  1978,  White,  1979),  the  particle  size 
parameter  is  given  only  e.xplicitly  ui  the  threshold  friction  velocity  which  varies  with  particle  size 
(Fig  4b)  All  three  of  these  equations  show  that  the  transport  rate  at  high  winds  becomes  independent 
of  size,  the  computed  results,  when  a  ranges  from  0.005  to  0.05  cm,  are  essentially  identical  for  u. 
exceeding  approximately  80  cm/s 

In  conclusion,  numerical  computations  show  that  most  of  the  sediment  transport  rate  equations  are  in 
reasonably  good  agreement  with  each  other  (within  a  factor  of  2  in  most  cases)  over  a  wide  range  of  wind 
velocities  and  sediment  sizes.  It  is  emphasized  that  such  close  agreement  among  various  formulas  is  quite 
remarkable  considering  the  fact  that  the  data  scatter  is  typically  a  few  orders  of  magnitude  in  most  reported 
experimental  data.  For  the  purpose  of  developing  an  algorithm  to  extract  the  size  distribution  information 
of  airborne  particles,  the  simple  expression  of  Bagnold  ( 1 94 1 )  is  selected. 

Although  there  are  several  studies  discussing  the  modification  of  transport  rate  due  to  particle  shape  factor, 
size  distribution,  and  sediment  density  (e  g.,  Williams,  1964,  Willetts  et  al.,  1982;  Willetts,  1983),  the 
proposed  "correction"  does  not  appear  to  be  conclusive  for  one  major  reason:  the  effects  of  those  param* 
eters  on  the  wind  friction  velocity  are  significant.  Much  of  the  deviation  be^vecn  the  transport  rate/wind 
friction  velocity  observations  and  the  basic  functional  form,  such  as  that  of  Bagnold's,  can  be  traced  to  the 
modification  of  u,.  This  will  be  further  explored  in  the  Phase  II  study. 
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4.  PARAMETERIZATION  TO  DERIVE  THE  SIZE  DISTRIBUTION 


To  investigate  the  feasibility  of  applying  the  sediment  transport  equations  to  denve  the  number  density  and 
size  distribution  of  particles  in  air,  an  equivalent  concentration,  C,  was  introduced.  Integration  over  the 
saltation  la>'er  of  the  product  of  C  and  the  particle  convective  velocity,  U^,  therefore  equals  the  sediment 
transport  rate.  Using  the  Bagnold  transport  rate  equation,  Eq.  (12)  w-as  generated  and  is  rewritten  here 


r  3 1  / 

=  nia,z)^  t/  (o.2)£far/r  =  1.81  — I  —uJ 

■''I  ”'‘’1  6  ^  V  j  5 


(39) 


where  is  the  median  diameter  of  the  bed  material.  The  square  bracketed  term  in  the  integrand  is  the 
equivalent  concentration.  To  allow  numerical  evaluation  of  the  vertical  size  distribution,  a  series  of 
parameterizations  is  required,  as  is  discussed  in  the  following. 


4.1  Size  Range  of  Mobilized  Sediment 

The  integration  limits  (a,.  Oj]  determine  the  size  range  of  mobilized  sediments  This  range  can  be  calcu¬ 
lated  using  the  threshold  critcnon  of  particle  motion  as  reported  by  Bagnold  (1 94 1,  see  also  discussion  of 
threshold  velocity  in  Section  3  I)  As  will  be  sho\^n  later,  the  effect  on  the  outcome  of  the  transport  rate 
computation  is  strongly  affected  by  but  not  o,.  The  lower  limit  of  the  size  integration  in  this  feasibility- 
study  is  set  to  zero 


4.2  Saltation  Layer  Thickness 

The  integration  range  (z,,  Zj)  represents  the  saltation  layer  thickness.  The  lower  limit  is  set  at  the  level  of 
one  sediment  diameter,  and  the  upper  limit  is  identically  the  saltation  height.  The  e.\act  magnitude  of  the 
saltation  height  is  not  well  deterrruned.  From  laboratory  results  (e  g,,  Bagnold,  1936,  Zingg,  1952,  White 
and  Schulz.  1977),  the  saltation  height  is  t>pically  reported  to  be  on  the  order  of  a  few  centimeters,  and  it 
increases  significantly  with  sediment  diameter  and  slightly  with  wind  velocity  The  unrealistically  low 
values  obtained  m  the  laboratory  are  due  to  the  geometric  limitation  of  the  expeiimenul  facilities.  In  the 
open  field,  such  as  described  in  the  reports  of  Bagnold  (1941,  p,  10)  cited  earlier  in  the  Introduction,  the 
saltation  height  is  on  the  order  of  approximately  2  meters,  it  is  larger  for  coarser  surfaces  and  smaller  when 
the  surface  is  composed  of  fine  sand.  Tiere  is  clearly  a  big  gap  between  the  field  conditions  and  laboratory- 
results  (ar.cJ  the  theories  developed  to  explain  the  laboratory  observations,  such  as  Owen,  1964,  WTute  and 
Schulz,  1977).  One  major  result  of  the  theoretical  studies  of  pa.nicle  trajectory-  shows  that  the  saltation 
height  is  closely  related  to  the  vertical  component  of  the  inii’al  liftoff  velocity,  w,.  The  liftoff  velocity-  can 
be  related  to  the  wjid  riiction  velocity.  In  Owen  (1964)  and  WTute  and  Schulz  (1977),  the  proportionality 
factor  derived  from  laboratory  observations  is  approximately  1.  Anderson  and  Haff  (1990)  suggested 
usii.g  »v,=5u.,  thus  Zj"  2!)u,-/2g,  which  is  closer  to  the  field  observation  value.  In  the  subsequent 
corip.jtations,  the  last  definition  for  Z2  is  used  to  set  the  upper  limit  of  integration.  It  is  noted,  however, 
that  tiie  choice  of  5  as  the  proportionality  constant  appears  to  be  arbitrary  and  that  the  resulting  expression 
show.s  no  dependence  of  saltation  height  on  sediment  diameter,  against  the  observational  evidence  described 
earlier.  Further  improvement  is  needed  in  the  future. 


4.3  Convective  Velocity 

The  convective  velocity  U^ia,  z)  is  assumed  to  be  the  wind  velocity.  L'(z): 

U(z)  =  2.5u.  In! 


V^o 


(40) 


where  the  dynamic  rouglincss  is  calculated  using  tiie  formula  z„=0  02065u.*/2jj  (Owen,  19o4,  .Anderson 
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and  HafT,  1990).  There  arc  few  reports  on  the  parameter  z),  but  the  numerical  compulation  of 
Anderson  (1986,  p.  27)  appears  to  support  that  is  proportional  to  V  over  an  extended  elevation. 


4.4  Vertical  Distribution  and  Shape  Factor 

The  function  n{a,  z)  is  expressed  as 

n(a,:)  =  Kn(a)/(z)  (41) 


where  is  a  scaling  factor  and  the  functions  n(a)  andyfz)  are  the  "shape  factors"  with  respect  to  size  and 
elevation.  From  extensive  research,  it  was  established  thaty^z)  attenuates  vertically  following  either  an 
exponential  law,  a  power  law,  or  a  combination  of  both.  In  either  case,  the  z>integration  can  be  carried  out 
explicitly.  Similarly,  there  are  a  few  choices  of  n(a).  The  two  most  frequently  adapted  distribution  forms 
are  the  log-normal  and  power  law.  Depending  on  the  fimctional  forms  specified,  the  scaling  factor  K  in  £q 
(41)  can  be  derived  from  Eq.  (39).  The  uitegrated  concentration  of  airborne  particles  within  a  defined 
volume  can  then  be  calculated  (see  Section  2  4).  To  illustrate  the  procedure  to  derive  size  distribution  from 
the  transport  rate  equation,  assume 


where  z,.  is  a  reference  elevation. 


n{a)  =  {alaj 


(42) 


(43) 


the  exponent  r=s  for  a/a„<\  (Region  I)  and  r-~l  for  a/a^>\  (Region  II),  where  both  s  and  /  are  positive 
real  numbers.  Observational  evidence  exists  showing  that  the  exponent  s  in  Region  I  (the  small-particle 
region)  is  somewhat  invariant  (e  g.,  Bagnold,  1941,  Figs.  41,  44)  at  various  stations  of  sampling  The 
exponent  /  in  the  large-particle  region  is  usually  larger  (shape  is  steeper;  I  igher  elevations  or  further 
downstream  of  the  sediment  source.  The  median  size,  a„,  of  the  airborne  particles  may  also  var\  with 
elevation  and  wind  condition.  Functional  expressions  for  the  variation  of  /  and  a„  with  z  or  u,  are  not 
available  at  this  stage,  but  data  from  Bagnold  (1941),  Chepil  and  Woodruff  (1957),  Williams  (1964),  and 
Nickling  (1975,  1978,  1983)  contain  reports  on  the  size  distribution  at  different  elevations  and  can  be  used 
for  establishing  the  relationship.  At  the  present,  the  parameters  1  and  a„  are  treated  as  constants. 


4.5  Analytical  Expression 


With  the  above  parameterization,  the  remaining  unknown  variable  is  the  scaling  factor  K  of  the  size 
distribu'  on  function.  By  choosing  a  seduiient  transport  rate  equation  (such  as  the  Bagnold's  transport  rate 
equation, ,  the  solution  for  K  cm  be  expressed  m  analytical  form 


where 


Q.72{am/aof^(u.^lg) 


V. 


m+  1 


Z2'"*'ln^-z,'""'li.  'i- 

Zq  Zq 


(44) 


(45) 


0.046 


u. 


2  'l 


4-1 


s  +  4 


go, 


-1 


m  y 


4-1 


(46) 


and  7(1=0  02065a, V2g,  z^=cl„,  and  Z2=25u.^l2g,  as  defined  earlier 
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5.  RESULTS  OF  SENSITIVITY  TESTS 


To  allow  numerical  evaluation  of  the  vertical  size  distnbution,  a  scries  of  parameterization  formulas  is 
required.  The  relevant  parameters  include  wind  friction  velocity  u.,  threshold  friction  velocity  u.p  fall 
velocity  size  range  of  mobilized  sediment  [a^,  a-^],  saltation  layer  thickness  (z,,  22],  convective  velocity 
Up(a,  2),  and  the  functional  representation  of  n(a,  2)  by  Kn{q)J{2).  For  the  purpose  of  demonstrating  the 
feasibility  of  obtaining  the  number  density  and  size  distribution  from  the  transport  rate  equations,  a  series 
of  analytical  expressions  was  established  for  the  various  parameters,  based  on  physical  considerations  and 
available  data  comparison.  Sensitivity  tests  were  performed  to  evaluate  and  rani:  the  factors  that  affect  the 
outcome  of  computed  number  density  of  airborne  particles.  A  set  of  computational  results  are  presented  in 
Figs.  7  through  9.  The  captions  provide  detailed  discussions  of  these  results.  Some  conclusions  are 
summarized  in  the  following: 

•  The  major  contribution  of  the  transport  is  from  the  large  particles  (Region  II,  This  is  because 

the  weight  of  transported  material  is  proportional  to  the  volume,  which  is  proportional  to  (see  panel 
(iv)  of  Figs.  7  through  9,  the  contribution  from  Region  11  is  denoted  with  0  and  the  total  contribution  is 
plotted  as  x], 

•  Because  the  transport  rate  is  dominantly  contributed  by  larger  particles,  the  vanation  with  the  large 
particle  distribution,  as  represented  by  the  parameter  /,  is  significant.  Notably,  for  a  narrow  size 
distribution  (/>3),  as  u,  increases,  the  size  range  of  mobilized  particles  and  the  number  density  increase 
(Fig.  7a(i)  and  (ii),  symbol  o).  However,  if  the  distribution  is  broad  (Fig.  7b)  while  the  size  range  still 
increases  with  increasing  wind,  the  number  density  may  actually  decrease.  This  is  due  to  the  fact  that 
the  motion  of  a  particle  derives  its  momentum  from  the  wind.  If  the  size  distribution  is  broad,  the 
number  of  larger  particles  mobilized  is  proportionally  increased  as  compared  to  the  case  of  a  narrow 
size  distribution.  Thus,  the  total  number  of  particles  needed  to  reach  equilibrium  with  wind  is  less  for  a 
broad  size  distribution.  Computation  shows  that  /=3  represents  a  critical  condition  (Fig  8a)  at  which 
the  total  number  of  particles  mobilized,  Y,  is  more  or  less  independent  of  the  wind  friction  velocity  (but 
the  representative  size  of  the  airborne  particles  increases  as  the  wind  velocity'  increases,  see  next  entry  ) 
For  narrower  distnbution  (/>3),  A'  increases  monotonically  with  u,,  for  broader  distribution  (/<3),  N 
decreases  as  u.  increases. 

•  Another  useful  indicator  of  the  size  distribution  of  airborne  material  is  the  representative  particle  size 
There  are  a  few  different  ways  to  represent  the  particle  size:  the  linear  dimension,  o,,  such  as  the 
diameter  of  sphencal  particles  or  the  axial  lengths  of  irregularly  shaped  particles,  the  square  root  of  the 
cross-sectional  area,  Oj,  or  the  cubic  root  of  the  volume,  03.  Each  expression  has  its  own  significance 
in  remote  sensing  applications.  For  example,  a,  is  relevant  in  measuring  the  transit  time  of  a  particle 
passuig  through  a  given  distance;  02  is  relevant  to  the  measurement  of  cross  sections  such  as  scattering 
and  visibility;  and  is  related  to  propertie.s  such  as  volume  absorption  and  mass  concentration.  As 
shown  in  panel  (iii)  of  Figs.  7  through  9,  the  characteristics  of  o,,  /=!,  2  and  3,  are  very  sensitive  to  the 
various  parameters  specifred,  especially  the  exponent  /.  Such  characteristics  can  be  used  to  develop 
methods  of  identify-ing  the  size  distribution  of  airborne  particles  by  combining  output  from  different 
sensors. 

•  The  effect  of  the  vertical  decay  function,  represented  by  the  parameter  m,  is  also  significant.  At  this 
stage,  it  is  found  that  the  computational  range  (in  terms  of  the  wind  friction  velocity  )  can  be  senously 
limited  for  a  particular  value  of  m.  For  example,  if  m=l.l,  no  computational  difficulties  were 
encountered  for  u,  up  to  200  cm/s;  when  /n=2,  it  was  found  that  the  computation  broke  down  near 
11. =80  cm/s  (see  Figs.  8a,  8b,  and  8c;  the  parameters  are  identical  to  those  of  Figs.  7a,  7b,  and  7c 
except  m)  The  coding  of  the  computational  programs  and  the  paramctcnzation  functions  can  be 
improved  to  e.xpand  the  range  of  applications. 
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Figure  7.  Case  studies  of  number  density  and  size  distribution  analysis  derived 

from  sediment  transport  rate  equations.  Plots  (a),  (b),  and  (c)  are  presented  to 
show  the  effect  of  parameter  I  (the  exponent  of  the  shape  function  of  large  particles) 
on  various  aspects  of  the  airborne  panicle  properties.  The  other  parameters  are 
ro=0.02065u.2/(2g);  22=(5u.)^/(2g);  o„=0.025;  j=4  1;  All  units  are  in  cgs. 

In  each  plot,  thiere  arc  four  panels,  (i)  through  (iv),  as  defined  below. 

(i)  The  nonnalized  particle  size  distribution  computed  with  u,  at  20,  40,  60,  200  cm/s.  The  increase 

of  the  size  range  of  mobilized  panicles  at  greater  u.  is  revealed  by  the  broadening  of  the  distribution 
curves.  Compare  panels  7a(i)  and  7b(i),  it  is  seen  that  in  the  former  case  (narrow  distribution  in  the 
source  matenal),  as  wind  increases,  the  airborne  panicles  increase  monotonically.  However,  if  the 
distribution  of  the  source  material  is  broad  as  in  the  latter  case,  although  the  size  range  of  particles 
mobilized  still  increases  monotonically,  it  takes  fewer  saltation  panicles  per  panicle  size-division  to 
"absorb"  tlie  wind  momentum  because  the  airborne  panicles  in  this  case  are  larger.  TTie  critical 
"broadness  inde.x"  is  /=3,  at  which  the  number  per  size-division  remains  approximately  invariant  with 
wind  [panel  7c(i)l. 

(ii)  The  scaling  coefficient,  K,  of  the  number  density  distribution  (x),  and  the  total  number,  N,  of  particles 
integrated  over  the  size  range  (o),  plotted  as  functions  of  wind  friction  velocity.  These  plots  give  a 
more  direct  apprehension  of  the  discussion  presented  in  (i).  For  the  two  lowest  wind  velocities,  the 
magnitudes  of  both  K  and  N  are  similar  in  the  three  panels  7a(ii),  7b(ii)  and  7c(ii)  because  the 
mobilized  particles  are  mainly  of  small  size,  the  distribution  function  in  the  small  particle  region  for 
the  three  computations  was  identical  (i.e.,  s=4.1).  As  wind  increases,  more  large  particles  are  sent 
into  the  air,  and  the  composition  of  airborne  panicles  is  affected  by  the  size  distribution  of  the  large 
particles  in  the  source  material.  In  general,  the  total  number  of  particles  is  expected  to  increase  with 
wind,  such  as  in  panel  7a(i).  However,  if  the  supply  of  larger  particles  is  abundant,  that  is,  broad 
distribution  in  the  source  function,  the  total  number  may  actually  decrease  with  increasmg  wind 
because  the  broadening  of  the  mobilized  size  range  introduces  more  large  particles  (as  compared  to  the 
case  of  narrow-distribution  source  matenal)  to  absorb  the  increased  wind  momentum.  Because  of  the 
presence  of  more  large  particles,  the  total  number  of  airborne  particles  required  to  reach  an 
equilibrium  with  wind  is  actually  reduced  Tfic  critical  broadness  index  is  /=3  [panel  7c(ii)],  at  which 
the  total  number  of  particles  airborne  is  maintained  approximately  constant. 

(iii)  As  a  result  of  the  difference  in  the  source  material  (i  .e.,  /),  the  composition  of  airborne  particles  can  be 
quite  different  under  identical  wuid  condition.  A  useful  index  to  denote  the  airborne  particle 
composition  is  the  representative  size,  which  can  be  the  linear  size,  computed  as  the  first  moment  of 
size  distribution  fiuiction  0|  and  plotted  as  the  lowest  curve  in  tliese  panels,  the  square  root  of  the 
cross-sectional  area,  computed  as  the  square  root  of  the  second  moment  and  plotted  as  the  middle 
curve,  or  the  cubic  root  of  the  volume,  computed  as  the  cubic  roo'  'f  the  third  moment  and  plotted 
as  the  top  curve.  For  narrow-distributed  source  material,  the  variation  of  o,,  ;=1,  2,  and  3,  at  different 
wind  conditions  is  small  [panel  7a(iii)],  for  broad  distribution,  the  representative  size  can  be 
significantly  different  at  different  wind  velocities  [panels  7b(iii)  and  7c(iii)]. 

(iv)  This  panel  shows  the  contribution  to  tlie  momentum  balance  from  the  large  particle  region  (o)  as 
compared  to  the  total  (x).  It  is  seen  that  except  at  the  very  low  wind  condition,  where  only  small 
particles  are  mobilized,  the  contribution  from  large  particles  dormnates  the  momentum  balance 
between  wind  and  sediment. 
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from  sediment  transport  rate  equations.  (Cont.) 
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Figure  8.  Case  studies  showing  the  effect  cf  the  vertical  decay  function,  repre¬ 
sented  by  the  exponent  m  of  the  function  f[z)=r'".  Plots  (a),  (b),  and  (c)  are 
computed  with  identical  parameters  as  those  shown  in  Fig.  7  except  that  m=2.  The 
plotted  properties  in  panels  (i),  (ii),  (iii),  and  (iv)  are  the  same  as  in  Fig.  7.  Because  of 
the  problem  of  a  singularity,  the  computation  range  of  u,  in  this  Sgure  is  20  to  80  cm/s 


Figure  8.  Case  studies  showing  the  effect  of  the  vertical  decay  function,  repre¬ 
sented  by  the  exponent  m  of  the  function  f(z)*r^.  (Cont.) 
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Figure  9.  Case  studies  showing  the  effect  of  the  median  diameter  of  the  source 


material,  as  represented  by  the  parameter  (a)  <7„-0.05  cm; 

(b)  .  ;„*0.0125  cm.  The  reference  condition  is  Fig.  7a,  a„=0. 025  cm.  The  most 
obvious  effect  due  to  this  parameter  is  in  the  rate  of  increase  of  particle  numbers 
as  a  function  of  wind  [compare  panels  9a(ii)  with  9b(ii)  and  7a(ii)].  The  smaller  the 
median  diameter,  the  faster  the  rate  of  increase. 
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The  effect  of  the  median  size  of  particles  on  the  number  density  is  shown  in  Figs.  7a,  9a,  and  9b,  in 
which,  s,  and  m  are  identical  and  only  is  varied  (0  025,  0.05,  and  0.0125  cm,  respectively).  The 
rate  of  increase  of  the  particle  number  with  the  wind  friction  velocity  is  found  to  be  reduced  as  a„ 
increases  (compare  panel  (ii)  of  the  three  figures).  This  result,  again,  is  a  consequence  of  the 
momentum  balance  governing  the  dynamics  of  saltation  and  wind  flow.  The  wind  momentum  is 
transferred  to  particles  through  the  drag  force,  which  is  proportional  to  a^,  the  mass  of  particles  is 
proportional  to  a^.  For  the  same  increment  of  wind  velocity,  the  increased  momentum  is  balanced  by  a 
given  mass  of  saltation  panicles,  which  can  be  from  a  large  number  of  small  panicles  or  a  small 
number  of  large  panicles.  Accordingly,  the  rate  of  increase  in  the  number  of  airborne  panicles 
decreases  as  the  panicle  size  increases. 
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6.  DISCUSSION 


6.1  Integrated  Number  Density  and  Angle  of  incidence 

In  the  lower  boundary  layer  near  the  ground,  the  vertical  gradient  of  airborne  particle  distribution  is  large. 
Therefore,  the  integrated  number  of  particles  inside  the  propagation  path  of  a  target  acquisition  system  is 
strongly  dependent  on  the  incident  angle.  Qualitatively,  it  is  expected  that  the  number  of  particles  at  a 
grazing  angle  will  be  considerably  larger  than  those  in  the  more  zenith-looking  configuratiotts.  The  spheri¬ 
cal  coordinate  system  is  used  to  obtain  the  integrated  number  of  particles  within  the  propagation  path. 
With  the  assumption  that  the  horizontal  distribution  of  airborne  particles  is  relatively  uniform  and  that  the 
vertical  distribution  can  be  expressed  asy(z)=r'",  the  integrated  number  of  particles  within  a  pencil  beam 
was  plotted  in  Fig.  2.  The  results  show  that  the  pa’ticle  number  increases  drastically  as  the  incident  angle 
increases  beyond  approximately  60**.  Depending  on  the  vertical  gradient  of  particle  distribution,  it  is  not 
uncommon  to  see  a  few  orders  of  magnitude  increase  from  a  10®  to  an  80°  incident  angle.  The  large  in¬ 
crease  of  airborne  panicles  in  the  propagation  path  can  cause  serious  deterioration  of  system  performance. 

6.2  Transport  Rate  Equations 

It  is  surprising  to  see  that  of  the  many  forms  of  sediment  transport  rate  equations,  the  calculated  results  are 
in  close  agreement  with  each  other,  in  most  cases  within  a  factor  of  2  (Fig  5).  This  is  truly  amazing 
considering  that  the  coefficients  of  those  equations  are  usually  derived  from  fitted  data  from  various 
sources  ranging  from  well-controlled  laboratory  experiments  to  complicated  field  conditions,  of  which, 
order  of  magrjtude  data  scatter  is  typical  These  equations  can  be  roughly  divided  into  two  major 
categories  (Fig.  6).  The  first  kind  shows  distinct  size  dependence,  such  as  the  equations  of  Bagnold  (1941), 
Zingg  (1952),  Owen  (1964),  and  Iversen  et  al.  (1976).  In  these  equations,  the  sediment  size  is  e.\plicitly 
incorporated  either  directly  as  a  or  as  the  fall  velocity.  In  the  second  category,  the  dependence  of  the 
calculated  transport  rate  on  sediment  size  is  implicitly  incorporated  in  the  threshold  velocity  parameter. 
Computational  results  indicate  that  the  effect  on  the  transport  rate  is  much  weaker,  especially  at  higher 
wind  velocities.  The  equations  of  Kawamura  (1951),  Lettau  and  Lctiau  (1978),  and  White  (1979)  fall  into 
the  latter  category.  Most  recent  equations  mcorporate  the  threshold  friction  velocity  in  the  formulation 
to  ensure  that  there  is  no  transport  rate  at  winds  less  than  the  threshold  magnitude.  The  equations  of 
Bagnold  (1941)  and  Zingg  (1952)  do  not  explicitly  specify  the  threshold  friction  velocity  in  the  formulation 
These  equations  are  understood  to  be  applicable  only  at  Although  .;quations  without  u.,  are  likely  to 

produce  larger  errors  at  lower  wind  velocities,  Bagnold's  equation  is  chosen  for  the  development  of  the 
algonthm  to  derive  size  distribution  from  the  transport  rate  equations,  due  to  its  simplicity  and  sound 
physical  basis  of  derivation.  The  resulting  algorithm,  however,  can  easily  be  modified  to  adapt  other 
transport  rate  equations. 

6.3  Size  Distribution  and  Sediment  Flu.x  Formuiation 

The  results  from  the  sensitivity  tests,  as  outlined  in  the  last  section,  demonstrated  that  the  number  density  of 
airborne  particles  can  be  derived  from  the  transport  rate  equations.  Because  of  the  large  number  of  factors 
involved  in  the  sediment  transport  process,  the  critical  step  to  an  accurate  prediction  of  the  particle  number 
density  is  the  parameterization  procedure.  As  expected,  some  of  the  factors  are  better  understood  than 
others.  For  example,  the  threshold  condition  is  more  or  less  established  and  the  size  range  of  particles  to  be 
mobilized  can  be  derived  as  a  consequence.  On  the  other  hand,  our  understanding  of  some  other  factors 
such  as  the  variation  of  Zj,  m,  i  and  with  respect  to  z  and  u,  is  only  qualitative.  A  more  quantitative 
relationship  can  only  be  derived  from  experimental  data  or  numerical  simulations  of  the  sediment  saltation 
process  A  small  set  of  experimental  data  and  numerical  results  has  been  collected  dunng  the  course  of  this 
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research.  These  include  the  laborator>’  data  of  Bagnold  (1936),  Zingg  (1952),  Williams  (1964),  Willetts  et 
al  (1982),  and  Willetts  (1983);  field  data  of  Hsu  (1971,  1973),  Gillette  and  Blifford  (1972),  Gillette  and 
Walker  (1977),  and  Nickling  (1975,  1978,  1983),  and  numerical  simulations  of  Anderson  (1986), 
Anderson  and  Hallet  (1986),  Anderson  and  Haff  (1988),  and  Werner  (1987,  1990).  The  present  param¬ 
eterization  scheme  can  be  improved  with  the  augmentation  of  the  above  data  sets.  However,  it  is  realized 
tiiat  air  the  data  sets  are  limited  to  particular  environmental  conditions,  and  continuous  expanding  of  the 
data  base  is  important.  Also,  impressive  progress  was  made  during  the  last  few  years  on  the  numerical 
simulation  of  eolian  transport  processes  by  Anderson,  Haff  and  their  co-workers  (sec  Haff  and  Anderson, 
1989,  and  Anderson  and  Haff,  1990,  for  a  more  detailed  discussion).  It  is  very  promising  that  with  a 
concerted  effort,  we  would  be  able  to  establish  a  reasonable  range  of  each  less-understood  parameter. 
From  such  information,  the  possible  range  of  integrated  particle  number  densities  within  the  propagation 
path  can  be  computed. 
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7.  RECOMMENDATIONS 

Based  on  this  feasibilit>'  study,  it  is  concluded  that  the  number  density  of  the  saltation  particle  size 
distribution  can  be  derived  from  sediment  transport  rate  formulas.  From  our  comparison  study,  it  is  found 
that  the  numerical  difference  among  formulas  of  various  complexities  is  small  compared  to  the  scatter  of 
experimental  data.  For  the  purpose  of  developing  the  algorithm  of  the  size  distribution  function,  simpler 
formulas  are  sufficient. 

Results  from  the  sensitivity  tests  indicate  that  the  number  density  of  saltation  particles  is  critically 
controlled  by  the  parameters  2^(0,  «.),  m{q,  u.),  I{z,  «.),  and  «,).  In  the  feasibility  study,  a  simple 
parameterization  scheme  was  used  to  reach  an  analytical  solution  The  parameterization  can  be  improved 
based  on  physical  considerations  and  experimental  data  (e  g.,  Bagnold,  1936,  1941;  Zingg,  1952;  Chepil, 
1957a,  1957b,  Williams,  1964;  Hsu,  1971,  1973,  Gillette  and  Blifford,  1972,  Gillette  arid  Walker,  1977, 
Nickling,  1975,  1978,  1983;  Horikavva  et  al.,  1983,  1984,  1986).  Reanalysis  of  these  data  sets  may 
provide  the  necessary  functional  relationship  and  the  likely  range  of  variation  of  each  parameter.  Thus  the 
upper  and  lower  limits  of  airborne  particles  can  be  calculated  for  output  to  the  application  programs  of 
target  acquisition  systems  to  enhance  their  performance.  Of  these  data  sets,  the  laboratory  results  of 
Williams  (1964)  and  the  field  data  of  Nickling  (1975,  1978,  1983)  provided  detailed  documentation  of  the 
size  distribution  as  functions  of  elevation  and  wind  friction  velocity  and  appear  to  be  most  promising  for 
our  purpose  of  establishing  the  parameterization  functions. 

Williams  (1964)  investigated  the  variation  of  size  distribution  of  sedimentary  particles  al  five  elevations 
(0.5,  2.5,  5.5,  9.5  and  15  5  cm)  in  a  wind  tunnel  9.14  m  long  with  a  35.5  x  35.5  cm  cross  section.  Hu-ee 
different  types  of  particles  (spheres,  sand  grains,  crushed  quartzite)  with  three  initial  size  distribuiions 
(uniform,  syinmetrical,  and  truncated)  were  uses  as  the  base  samples.  Experimental  runs  at  three  wind 
intensities  (approximate  u.  =  38*50,  78-92,  and  117-135  cm./s)  were  conducted  with  each  sample;  thus,  a 
combination  of  27  cases  were  obtained.  The  size  distribution  at  each  elevation  as  well  as  the  first  three 
moments,  representing  the  mean  diameter,  sorting  factor,  and  skewness,  were  reported  From  the  data  set, 
Williams  also  computed  the  vertical  distribution  of  sediment  concentration,  reported  to  decay  exponentially, 
and  transport  rate,  varying  as  a  power  of  u.,  with  the  exponent  being  2.7  for  crushed  quartzite,  3.4  for 
sand,  and  4. 1  for  spheres.  Detailed  results  including  the  size  distribution,  mean  diameter,  sorting  factor, 
and  skewness,  all  as  functions  of  vvind  and  elevation,  as  well  as  the  vertical  distribution  of  sediment 
concentration,  were  presented  in  graphic  form  in  the  onginal  paper.  The  data  set,  with  all  the  detailed 
measurements,  will  serve  as  our  fundamental  database  of  pa.''ametenzation.  However,  due  to  facility 
limitations,  the  vertical  scale  of  laboratory  data  sets  is  typically  1  to  2  orders  of  magnitude  smaller  than  the 
field  conditions  (for  example,  the  saltation  height  in  laboratory  experunents  is  typically  a  few  centimeters, 
while  in  the  field,  the  top  of  the  saltation  layer  is  1  to  2  m).  Field  data  and  numerical  simulation  may  help 
the  extrapolation  of  laboratory  data  to  field  conditions. 

The  field  results  of  Nickling  (1975,  1978,  1983)  include  data  of  creeping  (by  an  underground  trap), 
saltation  (collected  by  a  wind-vane  trap  with  an  orifice  1  m  high,  1  cm  wide,  and  the  lower  edge  0.05  cm 
above  surface),  suspension  (by  drawing  air  samples  through  filters),  and  wind  velocity  distribution  al  eight 
elevations  (0.5,  1,  1.5,  2,  4,  6,  9,  12  m).  The  ground  samples  were  collected  regularly  The  stati.stical 
properties  of  the  size  analysis  of  most  of  the  samples  were  also  reported.  Over  a  period  of  two  months, 
fifteen  storm  events  were  encountered.  Results  on  the  relative  coniribution  of  creeping,  saltation,  and 
suspension  to  the  total  transport  rate,  wind  dependence,  mean  size,  sorting  factor,  skewness  of  suspended 
samples  at  all  heights,  and  the  effect  of  moisture  and  salt  concentration  on  the  transport  rate  were 
documented  in  great  detail.  TTie  field  data  set  of  such  details  is  rather  unique  and  very  precious.  However, 
the  experimental  site  is  a  nver  della  with  very  limited  size  variability,  and  the  atmospheric  stability  condi¬ 
tion  appears  to  be  predominantly  unstable.  Because  of  the  fine  sand  size  (median  diameter  approximately 
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60  /fln),  the  contribution  of  the  susp'^nsion  load  to  the  total  transport  is  relatively  large,  more  than  60%  in  4 
of  the  15  storm  cases  (as  compared  to  20%  in  the  case  of  regular  dune  sand,  Bagnold,  1941). 

Not  to  be  overlooked  is  the  work  of  Bagnold,  summarized  in  his  classic  book  The  Physics  of  Blown  Sand, 
and  Desert  Dunes  (1941).  In  a  slightly  different  context,  Bagnold  investigated  the  size  di.stribution  of 
saltation  sediment,  including  the  mean  distribution,  the  temporal  variation  during  unsteady  wind  episodes, 
and  spatial  variation  over  nonuniform  surfaces.  It  was  found  that  the  size  distribution  at  any  stage  can  be 
expressed  as  a  power  law  n{a)={ala„y,  where  the  exponent  r  =  j  in  the  small  size  region  (Region  I,  a<aj, 
and  r  s  -y  in  the  large  size  region  (Region  II,  a>a^);  both  s  and  /  are  positive  numbers.  It  is  quite 
interesting  to  see  that  the  shape  function  in  Region  I  remains  relatively  constant  and  resembles  closely  that 
of  the  source  material.  This  property  is  important  for  our  development  of  the  transition  between  saltation 
to  suspension  at  the  top  of  the  saltation  layer.  The  shape  factor  in  Region  II  is  not  as  constant,  showing 
dependence  on  wind  and  space  or  time. 

In  order  to  establish  parameterization  to  cover  a  wide  range  of  source  material  parameters,  a  more  exten¬ 
sive  database  will  be  needed.  Tliere  have  been  many  data  sets  reported  in  the  sediment  literature  of 
different  degrees  of  environmental  parameter  measurements  and  different  levels  of  data  analysis  details 
Careful  classification  and  a  cataloging  system  need  to  be  established  at  an  early  stage  to  make  efficient  use 
of  these  precious  data  sets  In  any  event,  the  wntcr  is  confident  that  there  are  sufficient  research  results  for 
establishing  a  reasonable  range  for  each  parameter  of  interest.  For  the  established  range,  the  upper  and 
lower  limits  of  the  number  of  saltation  particles  under  given  wind  and  soil  conditions  can  be  defined. 
Program  output  can  also  include  information  such  as  the  most  likely  condition,  worst-  and  best-case 
scenanos,  confidence  interval  and  otlier  statistics  usefiil  to  the  system  operations  and  decision  making 

Due  to  the  time  limit,  only  the  theoretical  aspects  of  the  turbulent  suspension  process  were  revnewed 
Because  fine  dust  can  reach  a  much  higher  altitude,  the  effects  to  the  military  operations  can  be  significant 
Luckily,  our  undefstanding  of  the  suspension  process  is  much  more  satisfactory .  The  data  agreement  on 
the  vertical  distribution  with  vanous  theories  is  considerably  better  than  the  saltation  counterpart,  although 
the  absolute  magnitude  remains  a  difficult  quantity  to  p.irameterize.  Saltation,  occurring  at  the  lower 
region  of  the  boundary  layer,  serves  as  the  boundary  condition  of  suspension  From  continuity 
consideration,  the  reference  concentration  of  suspension  can  be  set  to  equal  the  saltation  "concentration"  at 
the  top  of  the  saltation  layer.  The  concept  of  matching  saltation  and  suspension  needs  to  be  tested  further 
with  expenmenDl  data. 

With  the  parametenzation  procedure  established,  programming  of  the  numencal  procedure  to  generate 
results  of  number  density  and  size  distribution  and  the  application  programs  designed  for  end  users  can 
proceed.  The  major  program  elements  and  application  considerations  arc: 

•  Develop  numerical  integration  subroutines  (quadrature)  such  tliat  the  computational  procedure  can 
handle  more  complicated  parameter  functions. 

•  Apply  the  saltation  results  to  the  suspension  component  of  the  sediment  transport  process.  Saltating 
particles  impact  the  surface  and  dislodge  small  dust  particles  that  can  be  suspended  into  high  elevations 
and  stay  aloft  for  a  long  period  of  time.  The  suspension  component  is  important  in  dutennining  the 
visibility  and  image  obscuration  at  high  elevations. 

•  Develop  the  coordinate  transformation  program  to  obtain  the  integrated  number  of  particles  in  the 
propagation  path  of  target  acquisition  systems.  Generate  results  in  graphic,  tabular,  and  computer  data 
file  formats  for  interfacing  with  other  remote  sensing  applications  or  additional  processing 

•  Develop  the  algorithms  that  relate  the  calculated  particle  size  concentration  to  visibility  ,  transmittance, 
absorption,  and  scattering  of  the  visible  to  infrared  optical  spectrum,  as  well  as  millimeter  radar  wave¬ 
lengths,  using  the  Mie  scaner  programs  available  at  QUEST  (counesy  of  Dr.  Yogesh  C.  Agrawal). 
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•  Generate  application  prograjns  on  a  PC  notebook  computer  to  generate  uitegrated  number  and  size 
spectra  of  particles  within  the  path  of  target  acquisition  systems  for  field  operation 

•  Compare  the  numcncai  output  witli  laboratory  and  field  data.  Design  and  carry  out  field  operation  and 
data  analysis  in  coordination  with  the  sponsor. 

•  Design  automatic  and  manual  modes  of  computation.  In  the  former,  the  program  requests  oidy 
minimum  uiput  parameters  such  as  estimates  of  wind  speed  and  soil  type.  The  program  will  generate 
results  based  on  a  set  of  prestored  parameters  to  provide  information  on  the  most  likely  conditions, 
worst-case  and  best<ase  scenarios,  and  assorted  statistics  to  assist  field  operations.  Manual  mode  is 
more  suitable  for  the  situation  when  more  detailed  environmental  measurements  are  available  and  is 
useful  for  planning,  forecasting,  and  hindcasting. 

•  The  appli>:.ition  programs  will  be  designed  to  be  operated  on  a  PC  notebook-based,  field-grade 
computer  User-friendliness  is  a  major  consideration  in  the  design  The  program  can  be  run  either  in 
Q&A  mode  or  it  can  be  menu  driven.  Results  will  be  presented  in  the  format  of  graphics,  tabulation, 
and  digital  files  for  communication  to  another  system. 

Undoubtedly,  this  list  of  elements  will  grow  as  the  research  progresses.  We  will  mainuam  continuous  com¬ 
munication  with  the  sponsor  to  update  and  improve  the  program  design  during  the  course  of  Phase  II  study. 
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